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Abstract 

The box ball system is studied in the crystal theory formulation. New con- 
served quantities and the phase shift of the soliton scattering are obtained 
by considering the energy function (or _ff-function) in the combinatorial 
i?-matrix. 



1 Introduction 

The box ball system (BBS) Q is one of the important classes of soliton 
cellular automata, namely integrable discrete dynamical systems with solitons. 
Recently, an interesting relation between the BBS and the crystal theory 
KMN| has been revealed. In particular, it is proved in TNS| that the scattering 



rule of two solitons in BBS is identical with the combinatorial i?-matrix |NY|. 



Motivated by this and the work |HIK|, the BBS is reformulated and generalized 
[HKT| in terms of the crystal theory. In this formulation, the time evolution 
is described by the row-to-row transfer matrix obtained as a product of the 
combinatorial i?-matrices. 

In the affine crystal theory, the combinatorial i?-matrix consists of the 'iso- 
morphism' and the 'energy function' or iJ-function. In the previous works, 
however, the _ff-function was not taken into account. In this paper, we study 
the BBS based on the crystal formulation and clarify the role of the iJ-function. 
The main results are as follows: 

• The phase shift of the soliton scattering is described by the i/-function. 

• New conserved quantities are obtained by including the i?-function into 
the row-to-row transfer matrices. 

In the simplest BBS, that corresponds to the n — 2 case in our paper, 
it is know n JITS I that a Young diagram obtained by the Schensted bumping 



procedure |F| is conserved under the time evolution. Their conserved quantities 
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are related to the number of solitons, and the quantities in this paper can 
be regarded as a generahzation to theirs. We note that equivalent conserved 
quantities are also obtained from those of the hungry Toda molecule equation 
through ultra-discretization |TN£],Q. 

The paper is organized as follows. In section 2, we recapitulate necessary 
facts from the crystal theory and reformulate the BBS. In section 3, we con- 
struct conserved quantities. The main theorem is given in section 4, where the 
scattering of solitons is studied. Section 5 is devoted for discussions. 



2 Preliminaries 

In this section we recapitulate necessary facts for the J7g(5[„)-crystal Bi and 
then reformulate the box ball system. 

2.1 Crystal Bi 

Fix an integer n G Z>2. For an arbitrary Z G Z>o we set 

Bi ^ {{vi,V2, - ■ ■ ,^i) I e {1,2, •••,"}, J^i < <••• < vi]- 

Bi can be identified with the set of semi-standard tableaux of shape [l) with let- 
ters in {1, 2, • • • , n}. For i = 0, — 1 we introduce operators e^, /i ; Bi — > 
Bi U {0} by the following rule. For i = 1, • • • , ?i — 1, ei{i/i, 1^2, ■ ■ ■ , t^i) is obtained 
by replacing the leftmost i + 1 in (vi, 1^2, ■ ■ ■ with i and fi{vi, 1^2, ■ ■ ■ ,1^1) by 
replacing the rightmost i with i + If we cannot find such letter to be replaced, 
the result of the action should be considered as 0. For i = the actions are 
given by 

eoii^i,i^2, ■ ■■ ,vi) ^ 5u^i(v2, • • • , vun), 
/o(i^i,J^2, • • ■ , i^i) = 5i^in(l, 1^1, • ■ ■ 

Notice if fih = b' for h, b' E Bi, then b = eib'. Bi is the crystal base Q of the 
Z-th symmetric tensor representation of the quantum affine algebra U^{sl„). For 
precise definitions, see ^, |KMN[ . 



Example 2.1 n = 4,1 — 5 

62(11334) = (11234) 
/2(11334) = 
/o(11334) - (11133) 

Drawing a line from b to b' with color i as & — ^ Bi turns out a colored 
oriented graph called crystal graph. 
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Example 2.2 (1) n: arbitrary 




(2) n=3 
B2 ■ 

11 — > 12 — > 22 

\ i \ i 
13 — > 23 

\ i 
33 

Here the colors of the arrows — >, i,\ are 1,2,0, respectively. 

2.2 Tensor product 

For b ^ Bi we define 

ei{b) = max{m > | e^^fe ^ 0}, ipi{b) = max{m > | ^ 0}. 

If 6 = {i/i, - ■ ■ e Bi, £i{b) (resp. 'Pi{b)) is tlic number of i + 1 (resp. i + nSio) 
in {yi, - ■ ■ For crystals Bi and Up the tensor product Bi (g) i?;/ = {6 (8) 6' | 
b & Bi,b' £ Bi'} is defined. The operators e^, /j act on Bii as 

e.(6 6') = ( (2.1) 

/.(6 6')=( ft^'; (2.2) 

[ if (^i(6) < £i(6'). ^ ^ 

Here b (E) and (8) &' are understood to be 0. Two crystals Bi i?;', -Be (X) 
are known to be isomorphic, i.e. the crystal graphs are the same. We denote 
the image of 6 (8> 6' under this isomorphism by b' (g) b. Namely, we have 

Bi ® Bi, ~ Bi, ® Bi (2.3) 
6 (8> 6' 6' (8> 6. 

We introduce so-called 'signature rule' to calculate the actions of e,,/j on 
the tensor product. It is most effective on the multi-component tensor product. 
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Consider an element &i (g) &2 C 
we associate an i-signature: 



)brn of Bl^ ®Bi. 



El . To this element 



.(0) _ „. 
. . (fc-1) 

pan' m 77- 



Set ri\ 



Vi{bi) ei{b2 



We construct rjl'^' from r]l ' by deleting an adjacent (H — ) 
(A; = l,2,---). Then we end up with the signature of form rj^ = 
+ (a, /? > 0). We call such signature reduced. Of course, there are 



many ways of deleting (H — ) pairs, but the reduced signature does not depend 
on them. Now the action of (resp. fi) on Bi^ ® Bi^ ■ ■ ■ (E) Bi^ is obtained 
by changing the rightmost — (resp. leftmost +) of ry, to + (resp. — ). If there 
is no sign to be changed, the action is understood as 0. It is easy to see that in 
the case of 2 components, this rule agrees with (^.l|),(2.2). For details, see next 
example. 



Example 2.3 n ~ 4:,m — 3 

Consider an element b = (1223) (112) (24) £ B4 
signature is given as follows. 



B^ Bo. The 1- 



b 

m 

The reduced signature is 



(1223) 



(112) 

- + + 



(24) 



1 1 



Vi = 



where the upper number signifies the component of the tensor product the sign 
belonged to. Therefore, we have 

eib = ei(1223) ® (112) ® (24)= (1123) ® (112) ® (24), 
fib = (1223) ® /i(112) ® (24)= (1223) ® (122) ® (24). 



2.3 Energy function 

Next consider a Z-valued function H on Bi® Bii satisfying the following prop- 
erty: For any b € Bi,b' ^ Bii and i such that ei{b eg) b') ^ 0, 



H{b ® 6') + 1 
i?(e.(6®6')) = { H{b®b')~\ 
H{b ® b') 



ifz = 0, Mb)>eoib'), Mb')>eoib), 
if i = 0, (po{b) < eo{b'), ipo{b') < eo{b), 
otherwise . 



H is known to exist and unique up to additive constant. The existences of the 
isomorphism (2.3) and energy function H are guaranteed by the existence of 
the i?- matrix. See [|KMN|. 
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Explicit rules to calculate the isomorphism and energy function are obtained 



in |NY]. Here we summarize the procedure to obtain them. We assume I > I'. 
Let b®b' be an element in Bi®Bii such as 6 = {vi, . . . ,vi) and h' = {y[, . . . , v[,). 
Let Xi (resp. yi) be the number of i in 6 (resp. b'). We represent h®h' hy the 
two column diagram. Each column has n rows, enumerated as 1 to n from the 
top to the bottom. We put Xi (resp. yi) letters i in the i-th row of the left (resp. 
right) column. The rule to obtain the energy function H and the isomorphism 
is as follows. 

(1) Pick any letter, say i, in the right column and connect it with a letter 
j in the left column by a line. The partner j is chosen so that j = 
max{i/fc I Vk < If there is no such j, we return to the bottom and the 
partner j is chosen so that j = max{;yfc}. In the latter case, we call such 
a pair or line "winding" . 



1 .1 



2 2 



unwinding 



winding 



(2) Repeat the procedure (1) for the remaining unconnected letters {I — 1)- 
times. 

(3) The isomorphism is obtained by sliding the remaining (k ~ I) unpaired 
letters in the left column to the right. 

(4) The value of the energy function is (-1) times the number of the "unwind- 
ing" pairs. 

Note that we normalized H so that we have H{{n'') (g) (n' )) = 0. Here and later, 
(n') means {n,n, ■ ■ ■ ,n) € Bi. When I' = 1, the rule simplifies to the following: 
If there exists k such that Vk < v' , then 



where j = max{fc | vy^ < v'}. Otherwise, 



{Vl, ■■■,Vl)® {v') ^ {Vl) ® {v', 1^1, • • • , Vl-l), 

H{{vi,---,vi)®{u'))^Q. 



(2.4) 



(2.5) 
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Example 2.4 1 = 4,1' ^2 



(1123) ® (23) ~ (12) 8) (1233) H = -2 
(1123) (12) ~ (13) ® (1122) H = -l 
(2344) (g) (12) ~ (44) (1223) H = 

2.4 Yang-Baxter equation 

Let us define the affinization AfF{Bi) of the crystal Bi. We introduce an inde- 
terminate z (the spectral parameter) and set 

AfF{Bi) = {z'^b\deZ,beBi}. 

Thus AfF{Bi) is an infinite set. z^b € AfF{Bi) will often be written as b. 

Definition 2.5 A combinatorial R-matrix for a crystal Bi S;/ is a map R : 
Aff{Bi) (g) Aff{Bi,) — > Aff{Bi,) (g) Aff{Bi) given by 

R{z% z'^'b') = 2'^'+"^''^''%' z'^-"^''^''%. 

The following result is a direct consequence of the ordinary [i.e. not combi- 
natorial) Yang-Baxter equation. 

Proposition 2.6 (Yang-Baxter equation) The following equation holds on 
Aff{Bi) (g) Aff{Bi,) (g> Aff{Bin). 

[R (g) 1)(1 (g) R){R (g> 1) = (1 i?)(i2 (g) 1)(1 i?) 

2.5 Box ball system 

Here and in what follows, we set B ^ Bi and consider the crystal _B®^ for 
sufficiently large L. The elements of B®^ we have in mind are of the following 
form. 

• • • (n) (n) • • ■ 0(n) (j^i) (^/g) • ■ • (i^k) (n) (n) • • • (n) • • • , 

where j^i, i^2j • • ■ ) ^fc G {!) 2, • • • , n}. Namely, relatively few elements are non (n), 
and almost all are (n). In the propositions below, we embed, if necessary, B^^ 
into B®^' {L < L') by 

6i • • • 6z, 1-^ 6i • • • 6i, (n) • • • (n) . 

^ V ' 
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Lemma 2.7 By iterating the isomorphism Bi®B ~ B®Bi, we consider a map 



Bi®B® 



iB 
'bL 



■B( 
61 



B®Bi 
)bL(g)b. 



Then there exists an integer Lq such that b — (n') for L > Lq. 

Taking sufficiently large L sucli that tlie above lemma liolds, we define a map 
Ti : B®^ — > B'^^ by 61 (X) ■ • ■ (g) 6l ^-> 61 (g) • ■ • ® 61,. 

Lemma 2.8 For a fixed element of B®^ , there exists an integer Iq such that 
Ti = Tig for any I >lo. 



Both lemmas are obvious from (2.4), (E 



Following HKT |, we reformulate the box ball system. We consider such 
an element of B'^^, oft en c alled a 'state', as described in the beginning of this 
subsection. Lemma 2.7 and 2.8 enables us to define an operator T = lim;_>oo Ti 
on the space of states. Application of T gives rise to a transition of state. Thus 
it can be regarded as describing a certain dynamical system, in which T plays 
the role of 'time evolution'. From the same reason, T; may also be viewed as 
a time evolution. (In this paper, time evolution means the one by T unless 
otherwise stated.) 

Under the condition assumed in Lemma 2.7, the sequence &i ® • • • Cg) 6^ is 
determined uniquely from hi ® ■■■ bL and b = (n'), since the isomorphisms 
Bi ® B ~ B ® Bi are bijective. Hence, the time evolutions Ti (1 > /) are 



invertible. The inverse map of the isomorphism ( p.4[ ),(2.5) is given explicitly by 
the following rule. 

If there exists k such that > v' , then 



(zy') (X) (z^i, ■ ■ ■ ,1^1) 1-^ {i^i, - ■ ■ ,v 
where j = min{fc | > v'}- Otherwise, 

{v') (g) {vi, ■■■ ,vi) ^ (z^2, • • • , 1^1, v') ® (vi). 



(2.6) 



(2.7) 



As we see in the next section, these time evolutions T; (Z > 1) form a 
commuting family of operators and play an important role in our paper. 



3 Conservation laws 

In this section we construct conserved quantities under the time evolutions. 
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3.1 Definition of Ei 

Fix sufficiently large L and consider a composition of the combinatorial R- 
matrices 

7^^ = i?L L+i • • ■ R23R12 ■■ Aff{Bi) ® AfF{Bf^ Aff{Bf'' ® Aff{Bi). 

Here Ra+i signifies that the i?-matrix acts on the i-th and (i + l)-th components 
of the tensor product. Applying TZi to an element {n}) ®p (p = 61 <8) ■ • ■ 
we have 

L 

where b'-j^ (0 < j < L) is defined by 

Bi® B ® ■ ■ ■ ® B ~ B®---®B ®Bi 

V ' ^ V ' 

3 j 
{n'-)®hi® ■ ■ ■ ®hj ^hi® ■ ■ ■ ®hj ® b^^'> . 

3.2 Conservation of Ei 

Proposition 3.1 As maps from AS {Bi)® AS {Bv)® AS {B)'^^ to AS{B)'^^ ® 
AS{Bi) ® AS[Bi), the following are identical. 

TZi'TZiRi2 ~ Rl+il+2T^iT^i' (3.1) 

Here Til (resp. Uv ) acts on AS {Bi) ® AS {B)'^^ (resp. AS{Bi,) ® AS{B)'>^^ ), 
and acts on the other component as identity. 

Proof. Use the Yang-Baxter equation repeatedly. I 

Theorem 3.2 For an element p e B®^, we have 

(1) TiTv{p)^TvTi{p). 

(2) Ei{Ti,{p)) = Ei{p). In particular, Ei{T{p)) = Ei{p). 

Proof. Consider an element (n') (g) (n' ) ® p {p — bi ® ■ ■ ■ ® b^) of AS{Bi) ® 
AS{Bi,) ® AS{B)®^. Apply both sides of (|]|) to it. Then one gets 

LHS = z"'+"4i ® • • • ® z"^+"4l z^'''-P\n^') z^^'^P^n^), 
RHS = z"'^+"4i ® • • • ® z"'^+"4l ® z^^''-P\n^') ® z^''^P\n^). 

Rerebi ® • • • ® 6l = Ti,Ti{p)M <S>---®bL^ TiTv{p) and p = Ti{p),p = Ti,{p). 
Hj, Hj, H'j, Hj (1 < j < L) are values of the energy functions, which are not 
necessary here. (|3.lD obviously ensures (1) and (2). I 
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4 Solitons 



4.1 Soliton 

A state of the following form is called an m-soliton state of length /i, Z2, • • ■ , ^m, 

...[^i] [h] ••• M (4.1) 

Here ..[I]., denotes a local configuration such as 

■ • ■ (8) n n (g) I'l ® ® ■ ■ ■ ® J^; '^n (g) n (g) ■ • ■ , (n > i/i > > • ■ ■ > > 1) , 
/ 

sandwiched by sufficiently many n's. (In this section, the parentheses in the 
state are omitted.) 

Lemma 4.1 Let p be a 1-soliton state of length I, then 

(1) The k-th conserved quantity of p is given by Ek{p) — min(fc,Z). 

(2) The state Tk{p) is obtained by the rightward shift by Ek{p) lattice steps. 

Proof. (1) Recall that the conserved quantity Ek is a sum of local i?- functions 
-H, = -i?(6(J-i) ® bj) e {0, 1}. By the rule (2.2) and (2.3), one has -Hj = 1 
if and only if j = J + 1, J + 2, . . . , J + min(fc, I), where J is the position of the 
last component vi in the tensor product B^. Hence E/^ = min(fc,^). Similarly, 
the statement (2) follows from the rule (2.2) and (2.3). I 

Example 4.2 The time evolutior^ of 1 soliton of length 3. Here and in what 
follows . means n (n> 3). 

....332 

332 

332 

332.... 

Definition 4.3 For any state p, the numbers Ni ~ Ni{p) (I — 1,2,...) are 
defined by 

Ei = J2^HkJ)Nk, Eo^O, 

k>l 

Ni = — + 2Ei — Ei^i. 



By Lemma 4.1, we have 

Proposition 4.4 For m-soliton state ((^), Ni is the number of solitons of 
length I, 

Ni = I Ij = I}. 

This proposition implies the stability of solitons, since the numbers Ei{p), and 
hence Ni{p), are conserved. 

^ The convention of the box ball game here is different from | TS , ^ in that, bigger numbers 
are transferred before smaller ones. 
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4.2 Scattering of solitons 

The following is an example of the time evolution (for < = 0, 1, . . . , 6) of a state 
which shows the scattering of three solitons of length 3,2 and 1. 

Example 4.5 

332. ..11. ..2 

. . .332. .11. .2 

33. .2112 

33. ..1221 

33. .1.. 221 

33.1. ...221... 

3.31 221 

We introduce a labeling of solitons of length I using Aff {Bi ) for the lower rank 
algebra C/^(s[„_i). (See also [ HKT | for such an identification.) Suppose there 



is a soliton of length I ■ ■ ■ i/i (g) V2 (S) ■ • ■ i^i (S) ■ ■ ■ {n > i/i > 1^2 > ■ ■ • > > I) 
at time t. Say it is at position 7(i), if (i^i) is in the 7(i)-th tensor component 



of B^^. From Lemma 4.1 (2), the position j{t) under the time evolution Tk 
is given by "f{t) = min(fc,^)< + 7 (min(fc,Z) is the velocity and 7 is the phase) 
unless it interacts with other solitons. To such a soliton we associate an element 
z-r (1^1,1^1.1, • ■ • , z^i) e AfFiBi) for U'^{5lr,^i). 

Now consider a state of m solitons illustrated as below. 

...[^i] [h] M 



We assume solitons are separated enough from each other and h > I2 > ■ ■ ■ > Im- 
Since longer solitons move faster, we can expect that the state turns out to be 



after sufficiently many time evolutions. (The proof of this fact is found in that 
of the main theorem.) We represent such a scattering process as 



z^'-^b'^® ■ ■ ■ ® z'^-^b'2® z^'^b'^ 



Here z'^^ bj , z"' are elements of AfF (Bi . ) under the identification in the previ- 



ous paragraph. With these notations, the scattering process in Example 4.5 is 
described as 

z°(2, 3, 3) (g) z-S(l, 1) ® z-"(2) ^ z-^i3) ® z-^{l,3) ® ^-^(l, 2, 2). (4.2) 

Let us recall some useful fact derived from representation theory. Note that 
Ulj{sln) contains C/q(s[„_i) as subalgebra. This fact can be translated into the 
language of crystals and guarantees that 

Tfe commutes with e^, (z = 1, 2, • ■ • , n - 2) on Aff(i?;J (g) ■ • ■ ® Aff(B/,,J. 
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Here the actions of ei, fi {i = 1, 2, • • • , n — 2) on the multi-component tensor 



product can be calculated using signature rule explained in section 2.2. By the 
actions, the power of z in an element of Aff{Bi^) is unaffected. We call this 
property J7g(s[„_i)-invariance. This property is also used to prove our theorem. 



For instance, if we admit (4.2), we have 



e2(z°(2, 3, 3) ® 1) ® z~"(2)) = z°(2, 2, 3) ® 1) ® z"^^(2) 

^ e2(z"*(3) ® z-^{l, 3) ® 2, 2)) = z-*(3) ® z-'^{l, 2) ® z-^{l, 2, 2). 

Let us explain the J7g(s[„_i)-invariance in terms of tableau combinatorics. 
Consider a state, read letters from right to left and remove all n's. Suppose 
the resulting sequence of letters is w = ii ■ ■ ■ (1 < ifc < n — 1). From such 
a sequence (or word) one can define a semi-standard tableau (• • ■ ((ii <— 12) 
13) • • ■ <— ig) by the 'row bumping' algorithm The C/g(5[„_i)-invariance 
assures that the resulting tableau is invariant under the time evolution. In the 



example (4.2), the words before and after the scattering are given by -u; = 211233 



and w' — 122133, and the corresponding tableaux are the same: 

112 3 3 
2 

Of course, this invariance is valid in the intermediate stage of scattering. 
The main results of this paper are the following. 

Theorem 4.6 (1) The two body scattering of solitons of length li and length 
h (h > h) under the time evolution Tr (r > I2) is described by the 
combinatorial R-matrix: 

AS{Bi,)®Aff{Bi,) 2i AS{Bi^)®AS{Bi,) 
z^'bi (S) z''^b2 ^ ^"^+'^62 <8) z'''-^, 

with 6 = 2l2 + H{bi ® 62)- 
(2) The scattering of solitons is factorized into two body scatterings. 

Some remarks may be in order. 

1) The combinatorial i?-matrix in Theorem 4.6 has an extra term 2^2 in the 
power of z. However, the Yang-Baxter equation (Proposition 2.6) holds as it is. 

2) Although we do not consider the cases where there exist solitons with same 
length, some part of the results can be generalized to such situations. Let us 
consider a state 

•■[fcMfc] ■•■••[2M2] [iMi] , 
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consisting of Nj solitons of length j {1 < j < k). By the same argument given in 
the following proof, we can show that the scattering of these solitons is factorized 
into scattering of two bunches {1 < j < i < k) 



We conjecture this scattering is described by the product of NiNj combinatorial 
R matrices AfF{Bi) >g)AfF{Bj) ~ AfF{Bj) Aff(Bi). This conjecture is trivially 
true if each solitons are separated enough. 

Example 4.7 Using this theorem, the scattering process is calculated as 

z"(2,3,3)(8)z"6(l,l)0z-"(2) ^ {z°(2,3,3)«)z-6(l,l)} ®z-"(2) 
^ {z-6+4(3^ 3) ^ ^o-Hi^ 1, 2)} z-"(2) - z-2(3, 3) ® {z-^{l, 1, 2) ® z""(2)} 
^ z-^{i, 3) ® {^""+'(1) ® z~^'\l, 2, 2)} = {z-2(3, 3) ® z-^°(l)} 2, 2) 

^ {z-i"+2(3) ® 3)1 ^ ^-5(^^ 2, 2) = z-«(3) 3) ^ ^-5(1, 2, 2). 

T/ie result is independent of the order of the factorization due to the Yang-Baxter 
equation. 



Proof of Theorem (1) Two soliton scattering rule: Due to the Uq{5ln-i)- 
invariance, it is enough to check the rule for the highest weight elements (1*) (8) 
(1', 2'^) € Bi® Bk+i {i > I + k), i.e. elements such that ej{bi ® ^2) — for any 
j — 1, • ■ • , n — 2. The corresponding state is given by 



Now consider the time evolution by T^+i+i. If j > I, (1*) moves with velocity 
k + I + 1 and (l', 2*^) with k + I. At some time, we arrive at the state 

••■(g)l(g)---®l(X)n(X)---®ri(8)2(X)---(X)2(g)l®---®l(g)---. 



By the rule (2.4)( p.5D , it is easy to see that after t time units (t < i ~ k — I) 



from this moment, the state becomes 
■ • ■ (g) nig) ■ ■ ■ n (g) 1 ■ - 1 ■ - (g) n 2 (g) ■ - (g) 2 1 ^ • • • . 

(fc+i + l)t i-t I k l+t 

After that {t > i — k — I), the solitons never interact again. 

We have investigated the time evolution by Tk+i+i. However, we have the 
equality 

^r^fc+i + l = Tk+l+l^r {r > k + I). 
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If a and b are sufHciently large, we can reduce the observation of the scattering by 
Tr in the right hand side to that by Tk+i+i in the left hand side, which we have 
just finished. It is immediate to check that (1') (g) (1^2'=) ^ (l'+'=) ® (^-'=,2'=) 
under the isomorphism Bi (g) ~ Bi+k ^ Bi. Let <5i (resp. 62) be the phase 
shift of the soHton of length i (resp. soliton of length {I + k)). They are given 

by 

51 = (k + I + - k - I) + {k + I) + I - {k + I + -k-l) = 2l + k, 

52 = {k + I + - k - I) - i - I - (k + -k-l) = -{21 + k). 

Since ^^((1*) (l', 2*^)) — -k, we have the desired result. 

(2) Factorization property: We illustrate the proof by a scattering of three 
solitons of length 3,2 and 1. Other cases are similar. Now the initial state is 

P=-[3]..[2] [1] 

We use a similar technique to (1). By applying the operator (a >> 0) we 
have scatterings of [2] [1] and [3] (g [1] . 

T2"(P) = [1] [3]..[2] 

Then by r| (6 >> 0) we have the third scattering of [3] ® [2]. 

Tl^T^^ip)^ [1] [2] [3] 

After this the time evolution T'^ acts trivially. On the other hand, these evolu- 
tion can be rewritten as 

In the right hand side, we first observe three soliton scattering caused by T'^ 
(c >> 0) and then T^T2 act trivially. Hence, the three soliton scattering in the 
right hand side is factorized into 3 x (two soliton scattering) in the left hand 
side. I 



5 Discussion 

In this paper we have studied the box ball system in the crystal base formulation 
and clarified the role of i/- functions. The function plays two important roles 
in the system: (1) to construct the conserved quantities and (2) to describe the 
phase shift of the scattering. Similar results are expected for the generalized 



models defined in |HKT| 



As far as we know, there is no direct proof of the equivalence between the 
original box ball systems and their crystal theoretic reformulation. This problem 
is under investigation. 

Another important feature of the box ball systems is their relation to non- 



linear integrable equations through the 'ultra discretization' |MSTTT, TTMS 
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rNS|. It is an interesting problem to understand the direct relation between 
such nonlinear integrable equations and crystal theory. 
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